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Abstract 

We report on a search for D^—D^ mixing made by studying the 'wrong-sign' 
process — > K~^7r~. The data come from an integrated luminosity of 
9.0 fb-i of e+e- colhsions at ^/s ^ lOGeV recorded with the CLEO II.V 
detector. We measure the time integrated rate of the 'wrong-sign' process 
K^ir^ relative to that of the Cabibbo-favored process — > K~^7r^ to 
be i? = (0.332 ± 0.040)%. We study D°^K+tt- as a function of decay 
time to distinguish direct doubly Cabibbo-suppressed decay from D^—D^ mix- 
ing. The amplitudes that describe D^ — D^ mixing, x' and y', are consistent 
with zero. At the 95% C.L. and without assumptions concerning charge-parity 
(CP) violating parameters, we find {l/2)x'^ < 0.041% and -5.8% <y' < 1.0%. 
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Studies of the evolution of a or B'^ into the respective anti-particle, a K'^ or B^, have 
guided the form and content of the Standard Model, and permitted useful estimates of the 
masses of the charm and top quark prior to their direct observation. In this Letter, 
we present the results of a search for the evolution of the into the Our principal 

motivation is to observe new physics outside the Standard Model. 

A D'^ can evolve into a through on-shell intermediate states, such as K^K~ with 
mass, rriK+K' ='m^D°, or through off-shell intermediate states, such as those that might be 
present due to new physics. We denote the amplitude through the former (latter) states by 
—iy (x), in units of T do/ 2 0. Many predictions for x in the amplitude have been 

made 0. The Standard Model contributions are suppressed to |x| ~ ta^v? 6c ~ 5% because 

decay is Cabibbo-favored; the GIM cancellation could further suppress |x| down to 
10~^ — 10~^. Many non-Standard Models predict > 1%. Contributions to x at this level 
could result from the presence of new particles with masses as high as 100-1000 TeV P,[lO . 
Signatures of new physics include |x| 3> ||/|, or charge-parity (CP) violating interference 
between x and y, or between x and a direct decay amplitude. In order to assess the origin 
of a — mixing signal, the effects described by y must be distinguished from those 
described by x. 

The wrong-sign (WS) process, — > i^T+vr", can proceed either through direct dou- 
bly Cabibbo-suppressed decay (DCSD), or through state- mixing followed by the Cabibbo- 
favored decay (CFD), ^ ^ K~^tt~. Both processes could contribute to the time 
integrated WS rate R = {f + /)/2, and the inclusive CP asymmetry A = (f — f) / (f + f), 
where / = r{D^ K^TT^)/r{D^ —>■ K^tt"), and / is defined by the application of charge- 
conjugation to /. 

To disentangle the processes that could contribute to ^ K^n^ , we study the distri- 
bution of WS final states as a function of the proper decay time t of the Z)°. We describe the 
proper decay time in units of the mean lifetime, Tdo = 415 ± 4fs The differential 
WS rate is [0,0 

r{t) ^ [Rd + 4R~Dy't + \{x'^+y'^)e]e-\ (1) 

The modified mixing amplitudes x' and y' in Eqn. |l| are given hj y' = y cos 6 — x sin S, 
x' = X cos 5+?/ sin 5 and Rm = ^ {x'^ + 2/^) = | {^'^ + v'^) where 5 is a possible strong phase 
between the DCSD and CFD amplitudes. The symbol Rd (Rm) represents the DCSD 



(mixing) rate relative to the CFD rate. There are theoretical arguments that 5 is small [jT4 
which have recently been questioned |]15 



The influence of each of x', y', and Rd on r{t) in Eqn. |I| is distinguishable. Such 
behavior is complementary to the time dependence of the decay rate to CP eigenstates such 
as ~*K^K~ that is primarily sensitive to or that of D*^ -^K^CrVp that is sensitive to 
Rm alone. 

We characterize the violation of CP in state-mixing, direct decay, and the interference 
between those two processes, respectively, by the real- valued parameters Am, Ad, and 0, 
where, to leading order, both x' and y' are scaled by (1 ± Am/2), Rd Rd (1 ±^d), 
5—^6±(f) in Eqn. ^ [|1^. The plus (minus) sign is used for an initial (D^)- 

Our data were accumulated between Feb. 1996 and Feb. 1999 from an integrated luminos- 
ity of 9.0 fb~^ of e^e~ collisions at a/s ~ 10 GeV provided by the Cornell Electron Storage 
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Ring (CESR). The data were taken with the CLEO II multipurpose detector |]T^, upgraded 
in 1995 when a sihcon vertex detector (SVX) was installed and the drift chamber gas 
was changed from argon-ethane to helium-propane. The upgraded configuration is named 
CLEO II.V. 

We reconstruct candidates for the decay sequence D*^ n^D^, K^n^. The 

charge of the slow pion (vr^ or tt~) identifies the charm state at t = as either or 
We require the D*~^ momentum, p^*, to exceed 2.2 GeV, and we require the to produce 
either the final state K~^tt~ (WS) or K~7i~^ (right-sign (RS)). The broad features of the 
reconstruction are similar to those employed in the recent CLEO measurement of the D 
meson lifetimes . 

The SVX provides precise measurement of the charged particle trajectories, or 'tracks,' 
in three dimensions |2^. We are thus able to refit the and vr" tracks with a requirement 
that they form a common vertex in three dimensions, and require that the confidence level 
(C.L.) of the refit exceed 0.01%. We use the trajectory of the K~^7i~ system and the position 
of the CESR luminous region to obtain the production point. We refit the vr^ track with 
a requirement that the trajectory intersect the production point, and require that the 
confidence level of the refit exceed 0.01%. 

We reconstruct the energy released in the D*^ ^n^D^ decay as Q = M*—M—m.,r, where 
M* is the reconstructed mass of the irfK^n' system, M is the reconstructed mass of the 
K^n^ system, and is the charged pion mass. The addition of the production point 
to the TT^ trajectory, as well as track-fitting improvements, yields the resolution aq = 190±6 
keV The use of helium-propane, in addition to improvements in track-fitting, yields 
the resolution ctm = 6.4 ± 0.1 MeV I^Tj. These resolutions are better than those of earlier 
and permit improved suppression of background processes. 



studies 22-25 



Candidates must pass two kinematic requirements designed to suppress backgrounds 
from TT+TT^, K~^K~, ^multibody, and from cross-feed between WS and 

RS decays. We evaluate the mass M for —>■ K^t{~ candidates under the three alternate 
hypotheses —>■ vr+vr", D'^ — »• K~^K~ , and — > 7r~^K~ . If any one of the three masses 
of the nominal mass 



falls within 4 a [26 



rri, the D'^—yK^Ti candidate is rejected. A 



conjugate requirement is made for the RS decays. The second kinematic requirement rejects 
asymmetric decays where the pion candidate has low momentum with the requirement 
that cos 6'* > —0.8 where 6* is the angle of the pion candidate in the rest frame with 
respect to the boost. The relative efficiency for the CFD to pass the two kinematic 
requirements is 84% and 91%, respectively. 

We require the specific ionization {dE/dx) measured in the drift chamber for each charged 
track agree to within So" of the expected value; this is a loose criterion, and we vary the 
dE/dx requirement for systematic studies. 

We reconstruct t using only the vertical {y) component of the flight distance of the 
candidate. This reconstruction is effective because the vertical extent of the e^e~ luminous 
region has ay = 7/im ||2^. The resolution on the decay point {xy,yy,z^) is typically 
40 /im in each dimension. We measure the centroid of the luminous region {xb, yt, z^) with 
hadronic events in blocks of data with integrated luminosities of several pb^^, and an error 
on yb that is less than 5/xm. We reconstruct t as t = M/py x — yi,) / {ctdo) where Py is the 
y component of the total momentum of the K~^7r~ system. The error in t, at, is typically 
0.4 (in lifetimes), although when the D° direction is near the horizontal plane, at can be 
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large; we reject 12% of the CFD by requiring at < 3/2. Studies of the plentiful RS sample 
allow us to determine our resolution function [|l9l, and show that biases are negligible for 
the WS results. 

Our signal for the WS process K^tt^ is shown in Fig. We determine the 

background levels by performing a fit to the two-dimensional region of < Q < lOMeV 
versus 1.76 <M< 1.97 GeV that has an area 135 times larger than our signal region. Event 
samples generated by the Monte Carlo (MC) method and fully simulated in our detector 
|28| corresponding to 90fb~^ of integrated luminosity are used to estimate the background 



shapes in the Q—M plane. The normalizations of the background components with distinct 
distributions in the Q—M plane are allowed to vary in the fit to the data. The background 
distributions and normalizations in the and samples are consistent and constrained 
to be identical. We describe the signal shape in the Q — M plane with the RS data that is 
within 4 a of the nominal CFD value. The results of the fit are displayed in Fig. ^ and 



summarized in Table I. 



TABLE I. Fitted event yields in a region of 2 a centered on the CFD Q and 
M values. The total number of candidates is 82. The estimated background is 
37.2 ± 1.8. The bottom row describes the normalization sample. 

Component Yield 



D^^K+n- (WS) 44.81^:? 

random 7r+ W-^K+t^- 16.0 ± 1.6 

e+e-^cc bkgd. 17.6 ± 0.8 

e~^e~ ^uu, dd, ss bkgd. 3.6 ± 0.4 

W^K+TT- (RS) 13527 ± 116 



The proper decay time distribution is shown in Fig. 2 for WS candidates that are within 
2cr 1^ of the CFD signal value in the Q—M plane. We performed maximum-likelihood fits 
in bins that are 1/20 of the lifetime. The background levels are constrained to the levels 
determined in the fit to the Q—M plane. We use the resolution function in t to describe the 
e~^e~ —>■ uu, dd, ss backgrounds, and an exponential, folded with the resolution function, to 
describe the e"'"e~ — >cc backgrounds. The distribution in t of the RS data is used to represent 
the random 7r+, —* K^n^ background [^. The WS signal is described by Eqn. |I|, either 
modified to describe all three forms of CP violation (Fit A), without modification to describe 
mixing alone (Fit B), or with the mixing parameters constrained to be zero (Fit C). The 
effect of our resolution is always included. 

The reliability of our fit depends upon the simulation of the decay time distribution of 
the background in the signal region. A comparison of the proper time for the data and MC 
samples for several sideband regions yields a = 4.4 for 8 degrees of freedom and supports 



the accuracy of the background simulation |21 



Our principal results concerning mixing are determined from Fit A. The one-dimensional, 
95% confidence intervals for x', y' and Rn, determined by an increase in negative log likeli- 
hood (— ln£) of 1.92, are given in Table II. The fits are consistent with an absence of both 
mixing and CP violation. The small change in likelihood when mixing and CP violation are 
allowed could be a statistical fluctuation, or an emerging signal. 
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TABLE II. Results of the fits to the — > K^tt decay time distribution. Fit A 
allows both D^-W mixing and CP violation. In Fit B, we constrain Am-, Ad, and 
to zero. In Fit C, we constrain x' and y' to zero, so R = Rd and A = Ad- The 
incremental change in — In £ for Fit B (Fit C) with respect to the Fit A (Fit B) is 
0.07 (1.57). 



Parameter 


Best Fit 95% C.L. 


Fit A 


Most General Fit 


Rd 


(0.48 ± 0.12 ± 0.04)% 0.24% < Rd < 0.71% 


v' 


(-2 5+}« ± 3)% -5 8% < V < 1 0% 


x' 


(0± 1.5 ±0.2)% \x'\<2.9% 


(l/2)x'2 


< 0.041% 




CP violating parameters 


Am 


0.23toil ± 0.01 No Limit 


Ad 


-Omtoir ± 0.01 -0.36 < Ad < 0.30 


sin (j) 


0.00 ±0.60 ±0.01 No Limit 


Fit B 


CP conserving fit 


Rd 


(0.47lg-i^ ± 0.04)% 0.24% < Rd < 0.69% 


y' 


{-2.3t\i ± 0.3)% -5.2% <y' < 0.2% 


x' 


(0± 1.5 ±0.2)% < 2.8% 


(l/2)s'2 


< 0.038% 


Fit C 


No-mixing fit R = Rd , A = Ad 


R 


(0.332l°:°i ± 0.040)% 


A 


-0.02i[J:^|^ ± 0.01 -0.43 <A< 0.34 


R/tan^ 9c 


(1.24l°:i±0.15) 
) (1.28l|]:i± 0.15 ± 0.03) X 10"^ 



We make contours in the two-dimensional plane of y' versus x' that contain the true 
value of x' and y' with 95% confidence, for Fit A and Fit B. The contour is where — ln£ 
increases from the best fit value by 3.0. All fit variables other than x' and y' are allowed to 
vary to best fit values at each point on the contour. The interior of the contour is the tightly 
cross-hatched region near the origin of Fig. |^. The limits are not substantially degraded 
when the most general CP violation is allowed, in part because our acceptance, unlike that 
of earlier studies ||2^j2^, is uniform as a function of decay time. 

Many classes of systematic error cancel due to the similarity of the events that comprise 
the numerators and denominators of R and A. The dominant systematic errors stem from 
potential misunderstanding of the shapes and acceptances for our backgrounds. We vary 
the selection criteria to estimate these systematic errors from the data. The level and com- 
position of the backgrounds are sensitive to the requirements on momentum magnitude and 
direction, and dE/dx of the charged particle trajectories and contribute ±0.018%, ±0.018% 
and ±0.026%, respectively, to the systematic error in Rd- We also include the statistical 



uncertainty on the MC determination of the proper time for the e~^e ^qq backgrounds ||29 
in the systematic error. We total systematic error on Rd, x' and y' of ±0.040%, 

±0.2% and ±0.3%, respectively. A study of detector-induced and event-reconstruction- 
induced asymmetries in CFD limits the relative systematic error on A to < 1%. 



7 



If we assume that S is small, then and we can indicate the impact of our work in 

limiting predictions of D^—D^ mixing from extensions to the Standard Model. The 95% C.L. 
interval for x from Fit A has some inconsistency with eighteen of the predictions tabulated 
in Ref. [0. A new model |^ invokes SUSY to account for the value of e'/e [0] and estimates 
mixing with x = 0.6%, somewhat below our sensitivity. Another analysis notes 
that SUSY could induce Ad^SO %, just below our sensitivity. 

In conclusion, our data are consistent with no D^—D^ mixing. We limit the mixing am- 
plitudes, x' and y', to be {l/2)x'^ < 0.041% and -5.8%<?/'<1.0% at the 95% C.L., without 
assumptions concerning CP violating parameters. We have observed 44.8lg7 candidates for 
the decay D^^K+n- corresponding to R = (0.332l|]:|j| ±0.040)%. We observe no evidence 
for CP violation. These results are a substantial advance in sensitivity to the phenomena 
that contribute to the wrong-sign process ^K^tt~ . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. We wish to acknowledge and thank the technical staff 
who contributed to the success of the CLEO II. V detector upgrade, including J. Cherwinka 
and J. Dobbins (Cornell); M. O'Neill (CRPP); M. Haney (Illinois); M. Studer and B. Wells 
(OSU); K. Arndt, D. Hale, and S. Kyre (UCSB). We thank Y. Nir for his help in the 
development of our description of CP violating effects. This work was supported by the 
National Science Foundation, the U.S. Department of Energy, Research Corporation, the 
Natural Sciences and Engineering Research Council of Canada, the A. P. Sloan Foundation, 
the Swiss National Science Foundation, and the Alexander von Humboldt Stiftung. 
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FIG. 1. Signal for the WS process K^tt . The data are the full circles with error bars, the 
projection of the fit for the signal is cross-hatched, and the projections of the fit for the backgrounds 



from charm and light quark production are single-hatched. For part a), M is within 2a |2f:] of the 
CFD value, and for b), Q is within 2a @ of the CFD value. 
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FIG. 2. Distribution in proper decay time for ^ K~^tt~ candidates, and the best fit of type 
A, described in Table II. Tlie data are sliown as tlie full circles with error bars. The cross-hatched 
region is the sum of the fit contribution from the direct ^ K~^it~ decay and the fit contribution 
from the destructive interference with mixing, which is shown in the region with single, vertical 
hatching. The fit contributions from backgrounds charm and light quark production are shown in 
single, diagonal hatching. 
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FIG. 3. Allowed regions, at 95% C.L., in the y' vs. x' plane. The entire kidney shaped region, 
filled with tight cross-hatching, is allowed under Fit A of Table II, while Fit B, in which CP violation 
is assumed, allows the smaller region, which is overlayed and filled with looser cross-hatching. The 
allowed regions from studies comparable to Fit A, using K~^K~ [^], for which we assume 

6 = 0, and ^ K~^l~Vi |31|, are shown as single-hatched regions. The Bayesian approach is used 

ni. 
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